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ABSTRACT 

We study the observational constraints on the cosmic evolution of the relationships between the 
massive black hole (MBH) mass (M,) and the stellar mass (M*.sph; or velocity dispersion a) of the 
host galaxy/spheroid. Assuming that the M, — M,.sph (or M, — a) relation evolves with redshift as 
cx (1 -\~ z)^ , the MBH mass density can be obtained from cither the observationally determined galaxy 
stellar mass functions or velocity dispersion distribution functions over redshift z ~ — 1.2 for any 
given r. The MBH mass density at different redshifts can also be inferred from the luminosity function 
of QSOs/AGNs provided known radiative efficiency e. By matching the MBH density inferred from 
galaxies to that obtained from QSOs/AGNs, we find that r = 0.64t°;^^ for the M, — M*^sph relation 

and F = -0.2lt[j;3| for the M, - a relation, and e = Cll^o pg. Our results suggest that the MBH 
mass growth precedes the bulge mass growth but the galaxy velocity dispersion does not increase with 
the mass growth of the bulge after the quench of nuclear activity, which is roughly consistent with the 
two-phase galaxy formation scenario proposed byiOscr ct al. (2 01^ in which a galaxy roughly double 
its masses after z = 1 due to accretion and minor mergers while its velocity dispersion drops slightly. 
Subject headings: black hole physics — galaxies: active — galaxies: evolution — quasars: general 



1. INTRODUCTION 

The masses of the massive black holes (MBH; 
M,) are tightly correlated with the properties of the 
spheroidal components of their host galaxies, such 

as th e velocity d ispersion a (e.g., iGiiltekin ct al.! 

200l iGebhardt et al 200ft iFerrarese fc M erritt 20 0ft 



Tremaine et al.l 120021: 

i 



minosity L. 



201 iD. t h e lu- 



[Graham et al.l 

e.g., iKormendv fc Richstond 119951 : 
20011) ■ a nd the stellar mass M^, 
1991 IMcLure fc DunlopI [2002 



IMcLure fc Dunlo. 

dMagorri an et alT _ 
[Marconi fc Hunt! 120031: [HaFing fc Rixl I2004D . which 
suggests a strong link between the growth of MBHs 
and the evolution of their host galaxies (or 
ularly the galaxy spheroids). Feedback 



partic- 
from the 

nuclear activities is proposed to be responsible for the 
establishment of these relationships, either through 
momentum or energy driven winds to s e lf-re gulate the 



Kina 2003: Wyithe fc Loebl 20031: 


Murrav et al.l [20051 


di Mattco et ai.l 20051: ICroton et al. 


[20061: [Bower et al.l 


20061: ISomerville et all I2008D. However, the detailed 



physics on how the feedback mechanism processes is still 
not clear. 

Observational determination of the cosmic evolu- 
tion of the relations between M, and M*_sph (or 
(T or ispii) may reveal important clues to the ori- 
gin of these rela tions and pu t constraints on feed 



IWoo et al.lj2008 : 



back mechanisrns ([Shields et al.[ 2003[: [Peng et al 



200 



Woo et al.1'2006': 'Trcu ct al."2007: 'Salviandcr ct al.''200 



Alexander ct al. 2008.; Somervillc 2009; 



Kisaka fc Koiim. 
20101; [Bennert et al 



Merloni et a l.l[201Ct [Sarria et al.! 
20ro n2011i; .Cisternas et al.l [2011f 



Schujze fc Wisotzkil[2011[ : [Portinari et al.[[2012l : [Li et alJ 
2012t) . A number of studies have shown that the 



evolves with redshift (or cosmic tirne), cx ( 1 -|- z)^ , where 
F ~ 0.68 - 2.1 (e.g., [Merloni et aH[20Tol; [Bennert et al] 
[20T1I) . There are also tentative observational evidences 
suggesting that th e M, — a relation ma y also evolve 
with redshift (e.g., [Woo et al.[ [20061 [2008[ ). In most of 
those studies, the MBH masses are derived by adopt- 
ing the virial mass estimators, which are based on the 
mass estimates of several dozen MBHs in nearby AGNs 
through the reverberation mapping technique and a cal- 
ibration of those masses to the M, — a relation ob- 
tained for nearby normal galaxies (e.g., [Onken et al.l 
[20041: [Graham et all [2011D . However, the MBH mass 
estimated from the virial mass estimators may suffer 
from some systematic biases due to v arious rea sons (e.g., 
Krolik 2001; Collin ct ah 2006: Nctzcr fc Marziani 2010; 
Kollatschnv fc Zetzl[ [2011t [Graham et al.l [20TT[ )rmore- 



M, — M^^sph relation in active galactic nuclei (AGNs) 



over, the MBHs in AGNs are still growing rapidly, unlike 
those in nearby quiescent galaxies. Therefore, it is not 
yet clear whether the cosmic evolution of the Af, — Af*^sph 
relation and the M, — tr relation found for AGNs are bi- 
ased or not and whether the relations for normal galaxies 
has a similar cosmic evolution as those for AGNs. 

In this paper, we adopt an alternative way to inves- 
tigate the cosmic evolution of the M, — M^^gph relation 
and the Af, — a relation in normal galaxies by matching 
the MBH mass density inferred from normal galaxies at 
different redshifts with that inferred from AGNs. The 
evolution of the M, — M*^sph and Af, — a relations for 
normal galaxies is assumed to follow a simple power-law 
form, i.e., cx (1 zY . In Section |21 we estimate the 
evolution of the MBH mass density in normal galaxies 
by the two ways: (i) using the stellar mass functions 
(SMFs) of galaxies determined by rece nt observations 
over t he redshift range 2: ~ — 1.2 (e.g.. [Bernardi et al.l 
[20101: [Ilbert et al.l [20101 ) . and (ii) using the velocity dis- 
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persion functions (VDFs) of galaxies at z ^ — 1.5 
(|Bernardi eFal] 120101: iBezanson at al.ll201lD . The evo- 
lution of the estimated MBH densitie s depends on the 
parameter F. According to the simple iSoltanI ()1982f ) ar- 
gument, the MBH mass density evolution can also be de- 
rived from the AGN luminosity functions (LPs), as shown 
in Section r2.3[ where the parameter F is not involved. By 
matching the MBH density evolution inferred from prop- 
erties of normal galaxies to that inferred from AGNs, the 
cosmic evolution of the M, — Aif*,sph (or M, — a) relation 
is then constrained in Section [3l Discussion and conclu- 
sions are given in Sections |4] and O 

Throughout the paper, we adopt the cosmological pa- 
rameters J?o = 70.5km s~\ 17a 0.726, and f^M = 0.274 
(jKomatsu et al.l [20091 ). The mass is in units of Mq and 
the velocity dispersion is in units of kms~^. Given a 
physical variable X (e.g., mass, logarithm of mass, veloc- 
ity distribution), the X function is denoted by rix(A', z) 
so that nj5s:(X, represents the comoving number 

density of the objects (e.g., MBHs or galaxies) with vari- 
able X in the range X — > AT -I- dX at redshift z. 

2. THE MASS DENSITY EVOLUTION OF MASSIVE BLACK 

HOLES 

2.1. The mass density of massive black holes inferred 
from stellar mass functions 

The mass of a MBH in the center of a nearby (z ~ 0) 
normal galaxy can be estimated through the M, —M^_sph 
relation at z = 0, i.e., 

(log M,) (M,,,ph; z = 0) = (8.20 ± 0.10) + (1.12 ± 0.06) 

x(logA4,,ph - 11), (1) 

where (log Af.) (M*^sph; z = 0) is the mean value 
of the logarithmic MBH masses for galaxies with 
spheroidal mass A'/H.,sph, and the intrinsic scatter of 
log M.(M<,. sph: z = ) around this mean value is 0.3de x 
(|Haring fc Rixl [200l see also iMcLure fc DunlopI [200l . 
Currently, there is still no consensus on the M, — MH.^sph 
relation for galaxies at high redshift. In general, the 
Af, — Af^.sph relation may evolve with redshift and the 
evolution may be simplified by 

(log Af.) (A//*,sph; z) = (log A-/.) (Af*,sph; ^ = 0) 

+riog(i-fz) (2) 

as assumed in a number of previou s works (e.g., 

iMerloni et "311120101 : iBennert et"aLll2011| ). where the pa- 
rameter r describes the significance of the evolution, and 
the intrinsic scatter of the relations is assumed not to 
evolve with redshift. 

The mass function of MBHs (BHMF) at redshift z may 
be estimated by adopting the A/, — A/*_sph relation (eq.[2|) 
and the SMF of spheroids at that redshift if the SMF can 
be observationally determined, i.e., 

ogA/.(logA/,,z) = y" niogA/.,,pi.(logAf*^sph,2:) X 

F(log A/.; (log A/.))d log A/,,sph(3) 

where niogM.,,ph (log A'/H.,sph, -z) is the SMF of spheroids 
at redshift z, P is the probability density function of 
logAf, around (logAf,) and is assumed to be normally 
distributed with a dispersion of 0.3 dex. The SMF of 
spheroids at redshift z can be estimated by using the 



SMFs for galaxies with different morphological types and 
the bulge-to-total mass ratios B/T, i.e., 



ni, 



ni, 



og A/.,.p„ (log A/H.,sph, z) = Yl "log Af.,,ph (log M^^spY,, z) 



/'(A^*,tot, z)niog Af.,tot (log Af*,tot, z) 



dlogAf^,,tot 
dlogAfH.,sph' 



(4) 

where Af*,tot = M^^sph/ {B / T), niogA/.,tpt(log Af*,tot, z) is 
the SMF for ah galaxies, nfogM. ^pi, (log A^*,sph, z) is the 
SMF for the spheroids of those galaxies with Hubble type 
i, /*(Af*_tot, z) is the fraction of galaxies with Hubble 
type i to all galaxies, and the summation is over galaxy 
morphological types from E, SO, Sa-Sb, Sc-Sd to Irr. 

The SMFs for galaxies at z ~ with different mor- 
phological types have be en obtained from the Sloan Dig- 
ital Sky Survey (SDSS; IBernardi et alllMol see Table 
B2 therein). And the bulge-to-tota l mass ratio B/T has 
been given bv lWeinzirl et all ()2009l see Table 1) for more 
than a hundred nearby galaxies with different Hubble 
types. In summary, the bulge-to-total mass ratios are 1, 
0.28±0.02, 0.46±0.05, 0.35±0.10, 0.22±0.08, 0.15±0.05^ 
and for E, SO, Sa, Sab, Sb, Sc-Sd, and Irr, respectivclyEI 
Adopting these observations, the SMF of spheroids and 
the BHMF at z = can be estimated. 

For galaxies at redshift z ~ 0.2 — 1.2, t he to - 
tal SMFs have been obtain ed bv lllbert et al.l (l 2010[) . 
Perez- Gonzalez et all ()2008l ). iFontana et a l. (2006), and 



Borch et al.l ()2006f ). These SMFs are usually obtained 



from deep surveys with small sky coverage and may suf- 
fer the cosmic variance. To avoid the cosmic variance, we 
adopt the average total SMFs ac cording to the SMF s es- 
timated in the above papers. In lllbert et al.l ()2010D , the 
SMFs of quiescent early type galaxies (E-I-SO) at differ- 
ent redshifts arc directly obtained (sec Table 2 therein); 
and wc obtain the SMFs of all the galaxies by summing 
the SMFs of early type galaxies and those of "interme- 
diate activity" and "high activity" galaxies (see the def- 
inition of "interm ediate activity" and "high activity" in 
lllbert et al.l 120101 and see their Table 3). We assume 
that the fraction of the early-type g alaxies at an y given 
mass is the same as that given by lllbert et all poTol ) 
and hereafter adopt this fraction, together with the av- 
erage total SMFs, to calculate the MBH mass function 
for early-ty pe galaxies. The LFs of galaxies were es- 
timated bv iZucca et al.l ()2006f ) for four different spec- 
tral types, which roughly correspond to the morpholog- 
ical types of E/SO, Sa-Sb, Sc-Sd, and Irr, respectively. 
For each type of galaxies, the mass-to-light ratio can 
be estimated through their average colors, for instance, 
log(Af*,tot/iB) = -0.942 -hl.737(B-y)-h 0.15 for early- 
type galaxies by adopting the Salpeter initial mass func- 
tion and log(A'/*,tot/iB) = -0.942 -I- 1.737(B-y)- 0.10 
for late-type galaxies by adopting the Chabericr initial 



^ IGraham fc WorlevI lj2008l ) obtain the ratio of K-band bulge 
luminosity/flux to disk luminosity /flux for a large number of disk 
galaxies in which the dust extinction is considered and corrected. 
Assuming a constant mass-to-light ratio, we also alt e rnativ ely 
adopt the B/D ratios inferred from [Graham fc WorlevI l|2008l 'l in 
our analysis, i.e., log(_B/D) = -0.54, -0.34, -0.54, -0.6, and -1.2 for 
SO, Sa, Sab, Sb, and Scd, respectively, and obtain the constraints 
(r,e) = (0.68 ±0.24, 0.10 ±0.01), which is c onsistent with that ob - 
tained by adopting the B/D ratios given bv lWeinzirl et al.l l|200l V 
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mass function (iBell fc De Jond [200lt iBell et alj [200l 
iBernardi et alll2010l: IChabriedl2003[ ) The B -V colors 
for different morpholo gical types of galaxies are given by 
iFukugita et alJ ()1995f ). The luminosity evolut ion can be 
correc ted for each type of galaxies according to lBell et alJ 
(j2003[) . T he LP for galaxies with different morphologi- 
cal types (jZucca et al.l 120061 ) can thus be converted to 
the SMFs. According to these SMFs, the relative abun- 
dance of different late-type galaxies (Sa-Sb, Sc-Sd, and 
Irr) over z ~ 0.2 — 1.2 can be obtained at any given 
M*,tot- The SMF of spheroids and the BHMF can then 
be estimated if the B/T \s averagely the same for galaxies 
with the same spectral type (and maybe correspondingly 
the same morphological type) but at different redshifts. 
Consequently, the mass density accreted onto MBHs with 
mass > A/, can be obtained by 

/•OO 

pT\z; > M.) = / (M:-Af.)^iogM;(logAf:,z)dlogAf:, 

J M, 

(5) 

if mergers of MBHs do not significantly contribute to the 
MBH growth and the see ds of th ose MBHs are smaller 
th an M, [see Eg. ( 29) iii lYu fc L u (2004) and Eq. (35) 
in lYu fc Tremaind (j2002| )] . And this mass density can 
be directly matched by that inferred from QSOs/AGNs 
(see Section [273]) . Hereafter, we set the lower limit of the 
MBH mass in the above integration to be IO^Mq unless 
otherwise stated, as the smallest mass for those nearby 
MBHs, of which the mass is well measured and adopted 
to determine the M, — M+^gph (or M, — a) relation, is 
~ IO^A/q. Given observationally well determined SMFs 
(or VDFs below), the MBH mass density estimated from 
the normal galaxies depends on the evolution parameter 
F in the M, — AfH.,spii (or M, — a) relation. 

2.2. The mass density of massive black holes inferred 
from velocity dispersion distribution functions 

The mass of a MBH can also be estimated through the 
M, — cr relation at redshift z — 0, i.e., 

(log A/.) (ct; z = 0) = (8.12 ± 0.08) + (4.24 ± 0.41) 

x(logcr-2.30), (6) 

with an intrinsic scatter of 0.44 dex (jGiiltekin et al.l 
|2009| ). Assuming an evolutionary form similar to equa- 
tion ([2]), i.e., 

(logAf.)(a;z) = (logAf.)(a;z = 0) 

+riog(l + z), (7) 

the BHMF at redshift z can be estimated through 



niog Ai. {log M,,z) 



(a,z)P(logA/.;(logA/.))da, 

(8) 

where nf^^{a,z) is the galaxy VDF at redshift z, P 
is the probability density function of logAf, around 
(log A^,) (cr; z) and is assumed to be normally distributed 
with a dispersion of 0.44 dex. T he VDF for local gala xies 
has been estimated from SDSS (|Bernardi et al.l[20Tnl see 
Table B4 therein). At higher redshift z ~ 0.3-1.5, 
the VDFs have been obtained from UKIDSS (United 
Kingdom Infrared Telescope Infrared Deep Sky Survey) 
Ultra-Deep Survey (UDS) and NEWFIRM (NOAO Ex- 
tremely Wide-Field Infrared Imager) Medium Band Sur- 
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Fig. 1. — The cosmic evolution of the comoving MBH mass den- 
sity. In the left panel, the red solid circles with errorbars represent 
the mass densities for those MBHs with mass > 10^ Mq, estimated 
from the SMFs of normal galaxies by assuming an evolution with 
F = 0.64 in the Af. — A/, gph relation (see cq. O, which are best 
matched by the solid line, while the red open circles represent the 
mass densities estimated by assuming a universal Af, — Aft^sph re- 
lation without evolution (i.e., F = 0). The solid line and the green 
shaded region represent the MBH mass density inferred from the 
hard X-ray AGN LF by assuming e = 0.11 and its Icr uncertainty. 
The right panel shows the best-fit parameter (F, e) (denoted by 
crosses) and the contours of the confidence levels for them, ob- 
tained by matching the MBH densities inferred from the normal 
galaxies to that inferred from AGNs. The contours enclose the 
68.3%, 95.4%, and 99.7% confidence regions on the joint distribu- 
tion of the two parameters. 



vey (NMBS) (|Bezanson et all I2011D . Adopting those 
VDFs, the BHMF and consequently the mass density 
accreted onto MBHs with mass larger than M, can also 
be estimated (see Eq. [5l). 

2.3. The mass density of massive black holes inferred 
from AGNs 

The MBH mass density at redshift z can also be in- 
ferred front the LF of AGNs according to the simple 
ISoltanI (|1982l ) argument as MBHs obtained their mass 
mainly through gas a ccretion ( e.g.. [ Salucci et ahl 119991: 
Yu fc T remai nd 120021: iMarcon i et al.l 120041: lYu fc Lul 



2004 120081: IShankar et al.ll2004 120091) . 



1 - e 



P. (2; > Af.) / dz / dLy / dG 



CyP{Cy\Ly)Ly4>{Ly,z) 



dt 
dz 



,(9) 



where Ly is the AGN F-band luminosity, (j){LY,z) is 
the AGN y-band luminosity function, e is the mass-to- 
energy conversion efficiency, and Cy = L^oi/Ly is the 
bolometric correction (BC) for the Y band, Lboi is the 
AGN bolometric luminosity, and P{CY\LY)dCY gives 
the probability of a BC to be in a range Cy — > Cy -l-dCy 
given a Ly . The hard X-ray LF of AGNs is adopted here 
because a significant number of obscured AGNs can be 
detected only in the hard X-ray band while missed in 
the optical surveys. Using the hard X-ray AGN LF and 
the corresponding BC, the MBH mass density can be 
estimated given a constant e. 

The hard X-ray AGN LF has been estimated 
by a number of authors based on the surveys by 
ASCA, Chandra, and XMM in the past decade (e.g. , 
Ueda et al. 2003; La Fr anca et al.ll200£ : ISilverman et al 



20081 lEbrero et al.i, 20091: lYencho et al J 120091 lAird et al 



In this paper, we adopt the late st 2-10 keV X-ray 
AGN LF obtained bv lAird et all (120101 . the LADE model 
in their Table 4) over the redshift range < z < 3.5 and 
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Fig. 2. — Legends are the same as those in Figure[T] except that 
the red solid (or red open) circles with errorbars are estimated from 
the VDFs of normal galaxies and by assuming an evolution with 
r = —0.21 (or no evolution with F = 0) in the M, — u relation, and 
the red solid circles arc best matched by the solid line estimated 
from the AGN LF by assuming e = 0.12. 

extrapolate it to higher rcdshiftQ Adopting other ver- 
sions of the 2-10 kcV AGN LF has httlc effects on the 
results presented in the following section. The BCs at 
the hard X-ray band (2-10 keV, denoted as Cx) have 
been found to be luminosity dependent and the BC at 
any giv en bolometric l umino sity has already b e en de - 
rived bv lMarconi et al.l (j2004[) and lHopkins et all (j2007D . 
To obtain the MBH mass density from the AGN LF, 
we need to estimat e the probabili t y dist ribution func- 
tion of P{Cy\Ly). iHopkins et all (|2007D obtained the 
probability distribution function of Cy as a function of 
L hoi- By adopting the same procedures as those done 
in IHopkins et al.l ( 2007| ). we fit the BCs by a log- normal 
distribution with the following parametert0, 



(logCx) = log 



Cl 



V lOiOL. 



fei 



C2 



lOiOL 



with (ci, fci, C2, ^2) 

scatter CTiogCx 
(0.24,0.06 



Ox 



0.01). Similar tolYu 



(10) 

(26.75,0.38,7.44,0.01), and the 
-I- (72 with {a\,^,a-2) = 
l3 (l2008l) . we then 



add a scatter of 0.15 dex t o include the X-ray yariabi l- 
ities in AGNs according to IVasudevan fc FabianI (j2007f ) . 
We adopt this fitting formula of Cx in the calculations of 
the MBH mass density. In the calculations of the MBH 
mass density from AGNs for MBHs with mass larger 
than Ajf, , a lower limit is needed to set to the luminosity. 

^ Since the Compton-thick AGNs are not included in the hard 
X-ray AGN LF, which may be a fraction of ~ 20% of the total AGN 
population fMalizia ct al .1120091) . the MBH mass densities inferred 
from the AGN LFs above may be underestimated by 20% and thus 
e for the best match is underestimated by a factor of ~ 1.2. 

^ We adopt the QSO spectral ener gy distribution (SED) model 
constructed by IHopkins et al.l II2007I ). For details, the template 
SED consists of a power law in optical-UV band, i.e., Li, oc v°' 
where Li, is the energy radiated per unit time per unit frequency 
at frequency v, with a = —0.44 for l/^m < A < ISOOAand 
a = -1.76 from 1200A to 500A. At the wavelength longer than 
A > Ifim, an infrared "bump" from reprocessing of optical-UV-X- 
ray emission is adopted and truncated as a Rayleigh- Jeans tail of 
the black-body emission {a = 2). The SED at energy high than 
0.5 keV is determined also by a power law, with a slope of —0.8, 
and an exponential cut-off at 500 keV. For any given monochro- 
matic luminosity at 2500A, the SED is renormalized to give the 
optical to X-ray slope of ctox = 0.384 log(L^{25ooA)/-f^i/{2keV)) 
and the points between 500 eV and 500A are connected with a 
power law. The value of Oox depends on luminosity as ciox = 
-0.107 log(L^(25ooA)/erg s"! Hz-i) + 1.739. After the SED is con- 
structed, we integrate the SED to obtain the 2 — 10 keV luminosity 
Lx and the bolometric luminosity, and hence the BCs. We fit t he 
BCs by a double power-law function as shown in Equation JTOj. 



As the typical Eddington ratio fo r low-lu minosity AGNs 
is ar ound 0.1 (jShen et al.ll2008D or 0.2 (iGraham et al.l 
1201 If ), therefore, it is reasonable to set the lower limit 
Lboi(M.) ~ 0.1Ledd(A/.) or 0.2Ledd(M.). Our calcu- 
lations show that the difference in the lower limit does 
not lead to significant difference in the results. 

3. CONSTRAINTS ON THE COSMIC EVOLUTION OF THE 
M. - M._sph RELATION AND THE M. - a RELATION 

In this Section, we obtain constraints on the evolution 
of the A/, — Af^.sph (or M, — a) relation by matching the 
MBH densities inferred from normal galaxies with that 
from AGNs using the standard statistics, i.e., 



E 



[pf^\z,,>M,)-pt^^{z,,>M.] 



(11) 



where Zi is the z-th redshift bin and the summation is over 
all the redshift bins, Spl''\zi,> A/.) and Spt™{z,,> 
M,) are the uncertainties in the MBH mass densities 
estimated from the normal galaxies and AGNs, respec- 
tively, considering the Icr errors in most of the fitting pa- 
rameters for the SMFs/VDFs and the AGN LFs, and the 
A/, — Af^.sph and M, — a relations. The errors in the nor- 
malizations of the A/, — AT^^sph and Af, — a relations and 
the AGN LF only introduce systematic shifts of all the 
estimates but do not lead to a change in the shape of the 
cosmic evolution of the MBH density, therefore, we do 
not include them in 5pT^{zi, > Af.) and 5p^'^^{z„ > Af.) 
in the fitting here. 

Figure [T] shows the mass densities of those MBHs with 
mass > lO^A^o estimated from the SMFs of normal 
galaxies, which are best matched by the MBH densi- 
ties estimated from the hard X-ray AGN LFs. The cor- 
responding parameters for the best match are (F, e) ~ 
(0.64lg'2g, O.IIIq'q^) and the errors are obtained for each 
parameter by marginalizing over the other parameter. 
According to our calculations, a F larger than 1.6 is 
excluded at the Str level. The value of F = 0.65 sug- 
gests that the AT, — Af^^gph relation evolves with red- 
shift positively, i.e., the relative positive offset of the 
A/, — A/*^sph relation at redshift z to that at local uni- 
verse increases with increasing redshift. Figure [5] shows 
the MBH mass densities estimated from the VDFs of 
normal galaxies, which are also best matched by the 
MBH densities estimated from the hard X-ray AGN LFs. 
The corresponding parameters for the best match are 
(F,e) = (-0.2lt°j8, 0.121° °^), which suggests that the 
Af, — a relation does not evolve with redshift. The e ob- 
tained by adopting the A/, — a relation is larger, which is 
because the MBH densities at ~ estimated from the 
A/, — cr relation is smaller than that from the Af, — Af^^sph 
relation. Considering the la errors in the normaliza- 
tions of the Af, — Af^.sph and Af, — a relations and the 
AGN LF, there are the additional errors ±0.03 in the 
estimated e. By combining these additional errors with 
the errors obtained from the above fitting and av- 
eraging the e obtained from the two relations, wc have 

I +0.04 

-o.mi 



obtained in lYu fc Lul ()2008[ ). 

As seen from Figures [T] and [2l the MBH mass den- 
sity at z ~ is about (3.5 — 4) x IO^Mq, which may 



be slightly smaller than that obtained by o thers (e.g., 
iGraham fc 671^120071 : iShankar et al.l [20091 ). This dif- 
ference is mainly due to the difference in the normaliza- 
tion of the adopted M, — <j relation and different treat- 
ment on the dependence of the MBH mass density on the 
Hubble constant. These differences lead to a change of 
the MBH density estimates at different redshifts by the 
same factor, but the shape of the MBH density evolu- 
tion does not change and thus the constraint on T is not 
affected. 

In the above calculations of the MBH mass densities, 
the SMFs and VDFs are extrapolated to the low-mass 
or low-velocity dispersion end, of which part may be 
not accurately determined by the obse rvations at all the 
redshifts considered i n this paper (see lllbert et all [2OIOI: 
iBezanson et~aI1l2011f ). To see whether the final results 
are affected by the extrapolation, we also set the lower 
limit of the MBH mass to IQ^Mq and re-do the above 

matching. And we find (r,e) = (0.61^^:2^,0.18^^:^2) to 
match the MBH mass densities estimated from the SMFs 
of normal galaxies with that estimated from the AGNs, 
and find (F, e) (-0.62l:°:20, 0.19lJ]:{]?) to match the 
MBH mass densities estimated from the VDFs of normal 
galaxies. The constraints on (F, e) for MBHs with mass 
> 1O^M0 are roughly consistent with that for MBHs 
with mass > IQ^Mq. 

4. DISCUSSION 

The positive evolution of the M, — Af^.sph relation 
found in this paper is consist ent with that found by 
Merloni et all (120101) (see also iJahnke et all 120091 and 
Bcnncrt e t al.l |2011[ ). which suggests that the growth 
of MBHs predates the assembly of spheroids and the 
spheroids experience additional growth after the quench 
of their central nuclear activities. The non-evolution of 
the M, — a relation found her e appears different from 
the positive evolution found bv lWoo et all (|2006l I2008D . 
We note here that this differen ce may be lessen ed as 
the MBH masse s estimated in I Woo et all (|2006f ) and 
IWoo et all (l2008h may be over-estimated by a factor of 
two as suggested by the uncertainties in the v irial fac- 
tor recently revealed by IGraham et al.l (j2011|) . How- 
ever, the exact reason for this difference is not clear 
as our results are obtained through the global evolu- 
tion of the MBH mass densities, different from the way 
adopted in iWoo et all 12006, 20081 ) for a sample of in- 
dividual AGNs / QSOs And the VDFs estimated by 
IBezanson et al.l (1201 ID are the very first estimates of 
VDFs at redshift z ^ 0, whi ch may suffer fr o m var ious 
uncertainties as discussed in IBezanson et ahl (j2011f ). If 
the non-evolution of the M^ — a relation is true, neverthe- 
less, it suggests the velocity dispersion of galactic bulges 
does not increase although the masses of the bulges can 
be significantly enlarged after the quench of their central 
nuclear activities. 

The evolution of the relations between the MBH mass 
and galaxy properties are intensively investigated theo- 
retically in the scenario of co-evolution of galaxies and 
MBHs since the discovery of these scaling relations. 
For those early models that adopt rapid and strong 
feedback due to energy output from the central AGNs 
which terminates star formation, the scaling relations 
are expected to evolve little and have very small scat- 
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ters fe.g.. [Granato et all 120041: iRobertson et all 120061: 
Idi Matteo et al.l 120051: iSpringel et al.l 120051 ). Later mod- 
els do suggest that the M, — M^^^sph relation evolve 
positively with redshift though with various degrees of 
evolution (e.g.. 'Hopkins et al. 2007; Croton et al. '20061 
Fontanot et al.|[2006; Malbon et al. 2007; Lamastra et all 
2010[ ). bv considering detailed dissipation processes oc- 
curred during major mergers of galaxies and the acquir- 
ing of bulge masses by dynamical processes such as disk 
instabilities or disrupting stellar discs. However, the ex- 
pected M» — a relatio n is almost independent of redshift 
(|Hopkins et al.l 120071 ). Apparently the constraints that 
we obtained in this pap er are roughly consi stent with 
the theoretical stu dies of iHopkins et al.l (120071). Further- 
more , we note that lOser et al.l (j2010l . see also lOser et al.l 
l2012f ) recently proposed a two-phase galaxy formation 
scenario, in which galaxies roughly double their masses 
after z ~ 1 due to accretion and minor mergers while 
velocity dispersion drops slightly. If MBHs mainly ob- 
tained their masses through efficient accretion triggered 
by major mergers, then our constraints are consistent 
with the two-phase galaxy formation scenario. 

4.1. Mass-to-energy conversion efficiency in 
AGNs/QSOs 

The constraints on F obtained above may depend on 
the use of a constant e, i.e., e is independent of M,, red- 
shift, and other physical quantities involved in the ac- 
cretion processes. We argue that the use of a constant e 
in AGNs/QSOs is appropriate as follows: (1) e is deter- 
mined mainly by the spins of MBHs in the standard disk 
accretion scenario; (2) the majority of the AGNs/QSOs 
are accreting via thin disks with high Eddington ratios 
(> 0.1); and (3) the spins of individual MBHs can quickly 
reach an equilibrium value and stay at that value for most 
of the AGN lifetime as suggested by theoretical models 
(e.g i Volonteri et III [200a ISImpiTa l20^ iHawlev et all 
120071 IMaio et al.ll2012D~ which suggests a roughly con- 
stant e at all redshift. Note that some authors introduced 
a dependence of e on t h e redshift or MBH mass, such 
as iDavis fc Laori (120111) : iMartfnez-Sansigre fc Rawlingsl 
(j2012f ): I Li et al.l (j2012l ). etc., according to the current 
observations. However, these res ults may be only due to 
some observational biases (e.g., iRaimundo et al.l |2012|) 
and need further investigation. Nevertheless, we note 
here that one could also introduce a cosmic evolution to e 
before understanding the underlying physics. For exam- 
ple, if we assume £(2) = max{0. 057, min[eo(l-l-z)'*, 0.31]} 
but the A/, — A'/*_sph relation does not evolve with red- 
shift, where 0.31 and 0.057 are the e that an efficiently 
accreting MBH-disk system could reach if the M BH spin 
is eith er 0.998 (the maximum spin of a MBH, see iThorni 
119741 ) or (a Schwarzschild MBH). In this case, an ac- 
ceptable fit can also be found and the best-fit parameters 
are (k, eo) = (1.3^°:^, 0.04^°:^^) for the M. -M,,,ph rela- 
tion, which suggests that e is significantly higher at 2: ~ 2 
than at z ~ 0. However, the MBH densities estimated 
from AGNs are an integration of over z, and is 
a function of e(z) and has to be determined to high red- 
shift. The MBH densities estimated from the SMFs in 
this paper only cover the redshift up to ^ 1.2 and may 
poorly constrain the evolution of e{z) at higher z. If e 
is significantly higher at higher redshift, a significantly 
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negative evolution in the M, — a relation is required. 
Future measurements on the SMFs and VDFs at z > 1.2 
may help to determine the MBH densities at higher z 
and thus may further help to put constraints on whether 
e significantly evolves with rcdshift. 

To close the discussion on this issue, we remark here 
that the constraint on the cosmic evolution of the M, — 
^jf*.sph (or M, — cr) relation obtained in this paper is ro- 
bust if the efficiency e of the efficient accretion processes 
in QSOs/AGNs is roughly a constant, which may be true 
as suggested by so me physical models of t he spin evolu- 
tion of MBHs (e.g IWonteri et al.ll2005t IShaDirol[200l 
iHawlev et al.ll2007l: [MaTo et al.ll2012ir 



4.2. Intrinsic scatters in the M, — Af,.sph or the M, — a 

relation 

The intrinsic scatters in the M, — A/*.sph and the M,—a 
relations have been assumed not to evolve with redshift 
in obtaining the constraints on the cosmic evolution of 
the relations. If the intrinsic scatter in the M, — M*^sph 
relation increases significantly with increasing redshift, 
the parameter F can still be consistent with to match 
the MBH densities estimated from normal galaxies to 
that from AGNs. To settle onto the observed local M, — 
-W*.sph relation with a smaller intrinsic scatter but with 
the same normalization, however, it is necessary for those 
galaxies at a fixed Af,^sph with relative large MBHs to 
accrete more stars and for those with relatively small 
MBHs not to accrete many stars after the quenching of 
the nuclear activities. This is not likely to be the case for 
the stochastic increasing of A/*^sph due to minor mergers 
or other dynamical processes like disk instabilities. 

In addition, we note that estimation of the intrinsic 
scatter of the M, — M*_sph (or M, — a) relation needs 
to determine the measurement errors in both the MBH 
and stellar masses, and it is still challenging to accurately 
determine these measurement errors (see discussions in 
[Graham et al. 2011). Graham (2012) shows that the to- 
tal scatter could range from 0.44 dex to 0.7 dex for dif- 
ferent types of galaxies (see Table 1 therein). 

4.3. Alternative M, — M^^gph or — a relation 

In the analysis in Section [H we adopt the single 
power-law form for the M, — Mj,,sDh re l ation given by 
iHaring fc Ril ([200l . Recently [GraE^ (|2012l) suggests 
that this relation may be better described by a broken 
power-law than the single power-law shown in Equa- 
tion (III). To see the effects of this new development 
on the constraints obtained above, here we replace the 
M, — M^,^sph relation at z = shown in Equation ([p 
by the broken power-law form given in iGrahamI ()20ia ). 
i.e., (logM.)(M,,sph;2 = 0) = (8.38 ± 0.17) + (1.92 ± 
0.38) log[M,.,ph/(7 X IO^OMq)] at A/*.,ph < 7 x IO^^Mq 
and (8.40±0.37) + (1.01 ±0.52) log[Af*,sph/(7xlOi"A//0)] 
at Af^^sph > 7 x IO^^Mq, respectively. Similarly, 
we also assume that the intrinsic scatter of this rela- 
tion is 0.3 dex. By doing the same analysis as that 
in Section [31 we obtain the constraints on (F, e) as 
{0A4:to^4,0mto^j). The best fit value of F = 0.44 is 
stih consistent with that obtained in Section[3|(0.64to;29) 
within 1 — a error but its uncertainty (^0 75) is large. 
Compared with the constraints on (F, e) obtained for the 



single power-law M, — Af*,sph relation, the larger uncer- 
tainty of F obtained here is mainly because the relative 
larger uncertainties in the slope of the broken power-law 
A/, — Af*,sph relation adopted here lead to larger un- 
certainties in the estimations of the MBH mass densities 
compared with that for the single power-law Af, — Af*,sph 
relation adopted above. The e obtained here (0.06) is 
also substantially smaller than that obtained for the sin- 
gle power-law Af, — Af^^sph relation, mainly because of 
the relat ively h igher zero-point of the broken power-law 
given bv lGraham (2012). 

In the analysis in Section [3l we also adopt the sin- 
gle power-law f orm f or the A/, — a relation given by 
iGiiltekin et al.l ()2009D . which is largely consistent with 
those e stimated by others ( see references therein). Re- 
cently, iGraham et al.l (j2011| ) updated the A/, — a rela- 
tion and found that this relation for the barred galaxies 
may be different from that for those non-barred galax- 
ies/ellipticals. If assuming that the Af, — tr relation 
is the same as that for ellipticals at cr > 180kms~^ 
[i.e., (logA/.)(cr;z = 0) = (8.22 ± 0.09) ± (5.30 ± 

0. 77) log(cr/200 km s^^) with an intrinsic scatter of 0.29] 
and the same as that for the barred galaxies at cr < 
180kms~i [i.e., (log A/.)(cr; z = 0) = (8.15 ± 0.06) ± 
(5.95 ± 0.44) log(cr/200 k ms~^) with an intrin sic scat- 
ter of 0.35; see Table 2 in IGraham et al.l 12011] . we find 
(F, e) = (-0.86tJJ jj, 0.14;°;°p. The constraint obtained 
here for F is substantially different from that obtained 
in Section [21 which suggests a significant negative evolu- 
tion of the velocity dispersion of individual big galaxies, 

1. e., the velocity dispersions of big galaxies decrease by 
~ 20%, which is marginally compatible with the hicrar- 
chical galax y for mation scena rio as recently proposed by 
lOser et al.l (e.g., I2010L l2012| ). However, the e obtained 
here is slightly higher than that obtained above for a 
single power-law Af, — a relation, which is mainly due 
to the smaller intrinsic scatter of the adopted Af, — cr 
relation and the smaller normalization for the relation at 
cr < 180 km s^^. Both of those factors lead to slightly 
smaller MBH mass densities in all the redshift bins and 
hence slightly higher e. 

5. CONCLUSIONS 

In this paper, we study the cosmic evolution of the 
A'jf, — A/,^sph relation and the Af, — a relation by a global 
method, independent of individual MBH mass estimates. 
We have estimated the cosmic evolution of MBH mass 
densities over the redshift range of z ~ — 1.2. The MBH 
mass densities are estimated from both the SMFs/ VDFs 
of normal galaxies determined by recent observations us- 
ing the A/, — A/,^sph and Af, — cr relat ions, and the 
AGN X-ray LFs according to the simple iSoltanI (|1982[ ) 
argument. By matching the MBH densities estimated 
from the normal galaxies with that from the AGN X-ray 
LFs, we obtain global constraints on the evolution of the 
A/, — A/,^sph relation and the M, — a relation. We find 
that the Af, — Af*_sph relation evolves with redshift pos- 
itively, i.e., oc (1 -h z)^ and F = 0.64to;29, though the 
significance level is not high; however, a F larger than 
1.6 is excluded at the 3cr level. We also find that the 
Ajf, — a relation appears not to positively evolve with 
rcdshift (F = -0.21^011). Our resuhs suggest that the 
MBH mass growth precedes the bulge mass growth but 
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the galaxy velocity dispersion does not increase with the 
mass growth of the bulge after the quench of nuclear 
activity, which is roughly consistent with the two-phase 
galaxy formation scenario proposed bv lOser et alJ ()2012l ) 
in which a galaxy roughly double its masses after z = 1 
due to accretion and minor mergers while its velocity 



dispersion drops slightly. 
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